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Abstract 

Glioblastoma multiforme (GBM) is a life-threatening 
brain tumor. Accumulating evidence suggests that eradi- 
cation of glioma stem-like cells (GSCs) in GBM is essen- 
tial to achieve cure. The transcription factor FOXMl has 
recently gained attention as a master regulator of mitotic 
progression of cancer cells in various organs. Here, we 
demonstrate that FOXMl forms a protein complex with 
the mitotic kinase MELK in GSCs, leading to phospho- 
rylation and activation of FOXMl in a MELK kinase-de- 
pendent manner. This MELK-dependent activation of 
FOXMl results in a subsequent increase in mitotic regu- 
latory genes in GSCs. MELK-driven FOXMl activation 
is regulated by the binding and subsequent trans-phos- 
phorylation of FOXMl by another kinase PLKl. Using 
mouse neural progenitor cells (NPCs), we found that 
transgenic expression of FOXMl enhances, while siRNA- 
mediated gene silencing diminishes neurosphere forma- 



tion, suggesting that FOXMl is required for NPC 
growth. During tumorigenesis, FOXMl expression 
sequentially increases as cells progress from NPCs, to 
pretumorigenic progenitors and GSCs. The antibiotic Sio- 
mycin A disrupts MELK-mediated FOXMl signaling 
with a greater sensitivity in GSC compared to neural 
stem cell. Treatment with the first-line chemotherapy 
agent for GBM, Temozolomide, paradoxically enriches 
for both FOXMl (+) and MELK (+) cells in GBM cells, 
and addition of Siomycin A to Temozolomide treatment 
in mice harboring GSC-derived intracranial tumors 
enhances the effects of the latter. Collectively, our data 
indicate that FOXMl signaling through its direct interac- 
tion with MELK regulates key mitotic genes in GSCs in 
a PLKl-dependent manner and thus, this protein com- 
plex is a potential therapeutic target for GBM. Stem 
Cells 2013:31:1051-1063 
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Introduction 



Glioblastoma multiforme (GBM) is the most common pri- 
mary malignant brain tumor and is highly aggressive and 
therapy-resistant [1-3]. Even patients with well-demarcated 
tumors in noneloquent areas that allow for gross-total resec- 
tion at surgery and who respond well to first-line chemo- 
therapies and radiotherapies frequently do not escape from 
subsequent recurrence and ultimately die from disease pro- 



gression. Therefore, there is an urgent need to develop 
novel therapies to effectively target resistant GBM cells. 
GBM is composed of mixed tumor cell populations includ- 
ing tumor cells with stem cell properties, termed glioma 
stem-like cells (GSCs) [4, 5] Accumulating evidence sug- 
gests that stem cell properties of GSCs contribute to thera- 
peutic resistance in GBM. Tumor stem cells are one, if not 
the only, cellular target in GBM, and therapeutic develop- 
ment targeting this subset of tumor cells may improve 
patient survival [6]. 
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The transcription factor, FOXMl is a member of the 
Forkhead box (FOX) family that consists of more than 50 
mammalian proteins with shared homology in the winged he- 
lix DNA binding domain. Many genes within this family play 
critical roles in cell cycle progression and cell fate decision. 
Initial studies of FoxMl characterized its role in liver devel- 
opment [7]. Intriguingly, FoxMl knockout mice fail to form 
hepatocellular carcinoma in a carcinogenic induction model, 
suggesting that FoxMl is necessary for tumor initiation in the 
liver [8]. Accumulated evidence suggests that FOXMl is a 
proto-oncogene with elevated expression in a number of 
human cancers such as liver, ovarian, breast, prostate, colon, 
and brain tumors including GBM. Downstream target path- 
ways of FOXMl include vascular endothelial growth factor 
(VEGF), matrix metalloproteinase-2 (MMP-2), and /^-catenin 
[9-11], each of which promotes tumor formation and progres- 
sion. FOXMl enhances cancer stem cell self -renewal through 
direct binding to /^-catenin inducing nuclear localization and 
transcriptional activity [10]. Taken together, abrogation of 
FOXMl signaling may provide multidirectional approaches 
for controlling cancers including GBM. 

Maternal embryonic leucine-zipper kinase (MELK) is a 
serine/threonine kinase and is abundantly expressed in GBM 
and various other cancers. Previously, we found that MELK 
is highly expressed in GSCs derived from GBM samples and 
its mRNA expression is inversely correlated with survival in 
GBM patients [12]. In addition, small inhibitory RNA 
(siRNA)-mediated MELK abrogation induces GSC apoptosis 
in vitro with less inhibitory effects on normal neural progeni- 
tor cells (NPCs) [12]. A recent study using a murine breast 
cancer initiation model indicated that ablation of Melk elimi- 
nates tumorigenesis in vivo; thus MELK is likely a therapeu- 
tic target not only for GBM but also for other cancers [13]. 

Prior studies have suggested that both FOXMl and 
MELK play essential roles in cell cycle progression, cancer 
cell growth, and maintenance of stem cell state of GBM. Fur- 
thermore, MELK and FOXMl are highly coexpressed and 
coregulated in GBM suggesting that they are functionally 
related. Therefore, we sought to test the hypothesis that 
MELK and FOXMl directly interact to regulate proliferation 
of GSC [14]. In this study, we interrogated the physical inter- 
action of FOXMl and MELK, as well as the functional roles 
of FOXMl-MELK signaling in GSCs. Here, we describe 
novel mechanistic evidence that FOXMl cooperates with 
MELK to regulate the mitotic transcriptome in GSCs. 



Materials and Methods 



Human Specimens and Tissue Culture 

Neurospheres (NS) derived from 10 GBM samples and 1 fetal 
brain were used in this study, as described previously [12, 15- 
17]. Regarding the nine samples collected at The Ohio State 
University (OSU), surgery was performed by I. Nakano and 
E.A. Chiocca in the Department of Neurological Surgery, and 
surgical specimens were processed for research under 
approved protocols (IRB Number 2005C0075). GBM157 was 
established in Dr. Kornblum's laboratory at UCLA, as 
described previously [18]. Established GSC cultures were cul- 
tured in defined medium containing Dulbecco's modified 
Eagle's medium (DMEM)/F12/Glutamax (Invitrogen, Carls- 
bad, CA, http://www.invitrogen.com) supplemented with 2% 
B27 supplements (Invitrogen) and 20 ng/ml epidermal growth 
factor (EGF) and fibroblast growth factor-2 (PeproTech, 
Rocky Hill, NJ, http://www.peprotech.com), as previously 



reported [19]. Mouse neural progenitors were cultured from 
the subventricular zone (SVZ) of the day 17 embryos, as 
described previously [17]. 293T cells were obtained from the 
American Type Culture Collection and cultured in DMEM 
(Invitrogen) supplemented with 10% fetal calf serum (Sigma 
Aldrich, St. Louis, MO, http://www.sigmaaldrich.com). For 
more detailed methods, see Supporting Information. 

Vectors 

All the FOXMl vectors (wild-type and mutants, S715A, 
S724A, double mutant S715/724A,S678A, T596A, TSAA,EE 
and the reporter 6x FoxMl) and PLKl (wild-type and domi- 
nant negative) were kindly provided by Dr. Zheng Fu of the 
Virginia Commonwealth University. All MELK vectors were 
described previously [20]. 

Phosphorylation Assay 

FOXMl cDNA was inserted in enhanced green fluorescent 
protein (pEGFP)-Cl (Clontech, Palo Alto, CA, http:// 
www.clontech.com) for mammalian cell expression and 
petl6b (Novagen) for bacterial expression. Flag-MELK and 
pet 16b-MELK 1-340 were described previously [20]. HEK293 
cells were transiently transfected with the indicated plasmids, 
and cell lysis was performed as described [20]. The com- 
plexes were trapped with GFP-trap beads (Chromotek, Ger- 
many, http://www.chromotek.com) and analyzed for the pres- 
ence of Flag-tag fused MELK with Flag antibody (Stratagene, 
La loUa, CA, http://www.stratagene.com). His-MELK 1-340, 
His-FOXMl, and His-cyclinA-CDK2 were expressed in bacte- 
ria and purified as described previously [20, 21]. His-FOXMl 
was phosphorylated in vitro by His-MELK 1-340, His-MELK 1- 
340D150A, or His-cyclinA-CDK2 for 1 hour at 30°C in a 
buffer containing 25 mM Tris at pH 7.5, 0.1 mM [}'-^¥] ATP, 
2 mM magnesium acetate, and 20 mM dithiothreitol (DTT). 

Gene Expression Omnibus Profile and TCGA Data 

Affymetrix Human Genome U133A Array was performed as 
described previously. The data have been submitted to the 
GEO database accession number CDS 18 15 [22, 23]. MELK 
and FOXMl expression (Affymetrix Human Genome U133A 
Array) data were downloaded from the CDS 1815 dataset and 
analyzed for grade III glioma and GBM. The Cancer Genome 
Atlas (TCGA) data are available through the TCGA Data Por- 
tal at http://tcga-data.nci.nih.gov [24]. 

Drug Treatment 

Siomycin A was obtained through the Developmental Thera- 
peutics Program NCI/NIH. Siomycin A and ON-01910 (PLKl 
inhibitor-selleck chemical Inc., Houston, TX, USA, http:// 
www.selleckchem.com) stock solution was prepared using 
dimethyl sulfoxide (DMSO) (Sigma Aldrich). In each experi- 
ment, DMSO alone was used in control samples at a concentra- 
tion in between 0.1% and 1% and identical to the concentrations 
of DMSO that was used in drug-treated cells in each experiment. 
We confirmed that the growth of the cells we used in this study 
is not significantly affected by DMSO at 1% or lower [18]. 

Statistics 

Quantitative data are presented as means ± SD, unless noted 
otherwise in the figure legend. The numbers of replicates are 
noted in the figure or legends. Comparison of mean values 
between multiple groups was evaluated by an ANOVA fol- 
lowed by Tukey's test. When multiple comparisons were 
involved. Holm's method was used to adjust the multiplicities 
to control the type I error rate that was less than 0.05 [25]. 
Spearman correlation coefficients (p) were used to assess the 
correlation between FOXMl and MELK expression in the 
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Figure 1. FoxMI expression is restricted in neural progenitor cells in the mouse brain. (A): Reverse transcription polymerase chain reaction 
(RT-PCR) analysis for FoxMI and Melk expression during brain development. (B): RT-PCR analysis for FoxMI and Melk in neural progenitor 
cells with or without differentiation. Expression is compared with the differentiation mai'kers for neurons (neurofilament heavy chain [Nfli]), 
astrocytes (Gfap), and oligodendrocytes (jvoteolipid protein [P/p]). For (A) and (B), Gapdh is used as internal control. (C): Left panels, top two 
images indicate immunohistochemistry for FoxMI (green) and Sox2 (red) in mouse brains of embryonic day 17 (E17). Lower four pictures, im- 
munofluorescence for FoxMI (green) and Tujl (red) in mouse brains of E17. Original magnification is xlO and x40. Scale bars = 100 and 20 
/(m. (a) CZ, (b) GZ, (c) ST represents magnified pictures of FOXMI (green) and Tuj-1 (red). Original magnification is xlO and x4D. Scale bars 
= 100 and 20 fim. Right panels represent immunohistochemistry for FoxMI (green), Gfap (red), BrdU (red), and Dcx in brains of the postnatal 
mouse day 30. Original magnification is x20-x40. Scale bars = 50—20 /(m, respectively. (D): Immunohistochemistry of FOXMI (brown) in 
adult mouse brain at the SVZ (left panel) and at the subgranular zone (right panel). Original magnification is x40. Scale bars = 20 /im. (E): Bar 
graph representing the relative number of NS f'omied from EGFP- or FoxMI -expressing NS derived from mouse E17 cerebral cortices. Asterisks 
(*) indicate statistical significance by t test. (F): Left panel exhibiting diagram of two different isofomis of FoxMI protein. Right panel indicating 
the relative number of NS formed from FoxMla or FoxMla/b siRNA transfected NS derived from mouse E17 cerebral cortices. NS forming 
assay was performed in triplicate in 96-well plate and repeated three times independently. Abbreviations: BrdU, Bromodeoxyuridine; CC, corpus 
callosum, CZ, cortical zone; DCX, doublecortin; DC, differentiated cells; EGFP, enhanced green fluorescent protein; ES, embryonic stem cells; 
E, embryonic brain; GFAP, glial fibrillary acidic protein; GZ, germinal zone; NS, neurospheres; P, postnatal brain; SVZ, subventrlcular zone; ST, 
striatum; V, ventricle; VZ, ventricular zone. 



GDS1815 microarray data and TCGA dataset [26]. Log-rank 
tests were used to compare the survival probabilities between 
groups in the xenograft mouse experiment. Comparison of 
mean values between two groups was evaluated by x" test or 
t test. Log-rank analysis is used to determine statistical signif- 
icance of Kaplan Meier survival curve. For all statistical 
methods, a p-value less than .05 was considered significant. 
For more detailed information, see Supporting Information. 



Results 



FoxMl Expression Is Restricted to NPCs in the 
Mouse Brain 

First, we investigated FoxMI expression during brain ontog- 
eny. Similar to Melk expression, FoxMI expression in the 
brain was predominantly detected during early and mid- 
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embryonic periods with a dramatic decline between embry- 
onic day 15 (E15) and E17 (Fig. lA). Expression of FoxMI 
in the adult brain was below detectable levels by reverse tran- 
scription polymerase chain reaction (RT-PCR). We then 
examined FoxMI expression in neural progenitor cultures 
grown as NS derived from embryonic mouse cortices. FoxMI 
expression was markedly high in proliferating NS, whereas its 
expression levels declined dramatically in the prodifferentia- 
tion — as indicated by expression of neuronal and glial 
markers — conditions within 6 hours, the earliest time point 
studied (Fig. IB). These data suggest that, similar to Melk 
expression, FoxMI is preferentially expressed in proliferating 
murine NPCs. 

We then performed immunohistochemistry of FoxMI with 
mouse embryonic and adult brains. FoxMI protein expression 
was restricted to the germinal zone (GZ) of embryonic day 17 
(E17) brains. The majority of FoxMI (+) cells were not colo- 
calized with neuronal marker TuJl and vice versa (Fig. IC, 
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left panels; Supporting Information Fig. SI). In contrast, we 
observed extensive overlap of FoxMl (+) cells with Sox2 
(+) cells in the GZ. In brains at the postnatal day 30 (P30), 
FoxMl expression was restricted to the proliferative zones 
lining the ventricles and was not colocalized with glial fibril- 
lary acidic protein (GFAP), suggesting that FoxMl is unlikely 
to be expressed in differentiated astrocytes or type B neural 
stem cells (NSCs) (Fig. IC, right panels). In addition, there 
was no detectable coexpression of FoxMl and the neuroblast 
marker Dcx. On the other hand, FoxMl was expressed by 
proliferating cells, as indicated by colocalization with BrdU. 
Taken together, these data suggest that the majority, if not all, 
FoxMl (+) cells are rapidly proliferating type C progenitor 
cells in the SVZ. When FoxMl expression was evaluated in 
the other neurogenic region, the subgranular zone (SGZ) of 
hippocampus, very few cells, if any, in the SGZ express 
FoxMl at P30 (Fig. ID). 

FoxMl Regulates NPC In Vitro 

To determine the function of FoxMl in NPCs, we assessed the 
effects of transgenic expression and silencing of the active 
FoxMl isoform b on NS derived from E17 cerebral cortices. 
These neural progenitors were transfected with either the expres- 
sion vectors or double-stranded RNAs designed to be siRNA for 
FoxMl (isoform b) or controls. For knockdown of FoxMl, we 
designed two control siRNAs that silence either nontarget 
sequence or the transcription repressor region of FoxMl isoform 
a, that is not present in the translationally active form of FoxMl, 
isoform b (Fig. IE) [27]. Overexpression of FoxMl (isoform b) 
resulted in a twofold enhancement of NS formation (Fig. IE; 
Supporting Information Fig. S2). Although siRNA-mediated 
silencing of the inactive FoxMl (isoform a) did not alter the 
growth of NS compared with the nontarget control, siRNA 
directed against the common region of FoxMl (isoforms a and 
b) resulted in a significant decrease in NS numbers (Fig. IF, right 
panel; Supporting Information Fig. S2). Collectively, these data 
suggest that FoxMl is essential for the growth of murine NPCs 
derived from embryonic brains in vitro. 

FoxMl Expression Is Markedly Elevated During 
Gliomagenesis 

Accumulating evidence suggests that GSCs and their somatic 
counterparts, NPCs, share various fundamental signaling path- 
ways to regulate their self-renewal and proliferation. We 
therefore sought to determine the role of FoxMl in GBM-like 
tumors and GSCs in mice harboring conditional inactivation 
of the human GBM-relevant tumor suppressors p53, Nfl, and 
Pten [28]. In this model, Nestin promoter-driven Bromodeox- 
yuridine (BrdU)(+) cells along the ventricular wall are the 
cells of origin for the cortical GBM-like tumors formed in 
adult mice [28]. When we compared FoxMl staining of the 
SVZ of wild-type mice and genetically engineered mutant 
GBM mouse brains, we observed sequential elevation of 
FoxMl expression along the ventricular wall during the 
course of tumorigenesis (Fig. 2A; Supporting Information Fig. 
S3). When FoxMl expression was compared between normal 
brain regions versus tumor areas of individual mice, tumors 
had more than 200-fold higher expression than normal brain 
regions (Fig. 2B). We then compared FoxMl expression in 
NS derived from the SVZ of wild-type mice and GBM-like 
tumors. FoxMl mRNA levels were 2.5-fold greater in GBM 
NS than in normal NS, indicating that FoxMl expression is 
elevated in GSCs compared to NPCs in vitro (Fig. 2C) [29]. 

The FOX family transcription factors play a crucial role 
in organ development and cancer initiation and propagation in 
a context-dependent manner [30, 31]. To determine which 
members of the FOX gene family are upregulated in GSCs, 



we evaluated expression levels of all the FOX genes in the 
transcriptome microarray data derived from 30 patient-derived 
GSC samples as well as 4 normal astrocytes primary cultures 
(Fig. 2D). Among 51 FOX family members, FOXMl exhib- 
ited strikingly higher expression in GSCs in comparison to 
normal astrocytes, unlike any other FOX family members. In 
all of the tested 10 GSC samples, we observed 6-12-fold 
higher FOXMl expression by quantitative RT-PCR compared 
to astrocytes (Fig. 2E). CD133 expression on the cell surface 
is associated with cancers stem cells in some, if not all, GBM 
tumors [1,4, 16, 32-34]. We further performed cell sorting by 
using CD133 marker and separated CD133 (+) and CD133 
(-) cells for GBM157 and GBM528 GSCs. qRT-PCR demon- 
strated higher expression of FOXMl and MELK in CD133 
(+) cells than CD133 (-) cells in both samples (Fig. 2F). 

We then examined FOXMl protein expression in human 
GSCs. In agreement with previous studies, immunohistochemis- 
try of a human GBM specimen demonstrated strong immunore- 
activity for FOXMl in the nuclei of GBM cells (Supporting In- 
formation Fig. S4) [35]. Immunocytochemistry of dissociated 
GBM NS demonstrated that FOXMl was colocalized with the 
stem cell-associated mai'kers including Nestin and SOX2 but not 
with the differentiation markers, TuJl or GFAP (Fig. 2G). Pro- 
tein expression analysis by flow cytometry also yielded similar 
results. Similar to the higher proportion of tumor cells expressing 
the NSC-associated protein SOX2 and in contrast to the lower 
proportion of tumor cells expressing the astrocytic marker GFAP 
in GSCs in serum-free conditions, FOXMl (+) cells were pre- 
dominantly detected in GSCs in serum-free medium (GBM 
NSs), when compared to GBM cells propagated in serum-con- 
taining medium (SPGCs) (Fig. 2H; Supporting Information Fig. 
S5). Subsequently, we analyzed FOXMl promoter activity in 
GSCs and non-GSCs in GBM30. GSCs manifested a signifi- 
cantly higher FOXMl promoter activity in comparison to SPGCs 
(Supporting Information Fig. S6). Taken together, FOXMl 
expression and activity are substantially elevated in GSCs com- 
pared to SPGCs and normal cells in mouse and human brains. 

FOXMl Is a Substrate for MELK 

Prior studies have suggested that both FOXMl and MELK 
play essential roles in the cell cycle progression, cancer cell 
growth, and maintenance of stem cell state of GBM [18, 20, 
32, 35, 36]. Nonetheless, no studies have demonstrated a 
direct molecular interaction of MELK and FOXMl. We there- 
fore sought to elucidate the potential molecular interaction of 
FOXMl and MFLK. First, we used the TCGA array data set 
(GEO database accession number CDS 18 15) to determine the 
similarities of the expression profiles between FOXMl and 
MELK in high-grade glioma (Fig. 3A). Using 24 grade III gli- 
oma samples and 56 GBM samples, we found that FOXMl 
RNA expression in individual tumors was correlated with 
MELK expression (p<.001). Significant correlation was also 
observed in newly diagnosed and recurrent GBM patient sam- 
ples at RNA levels (Supporting Information Fig. S7). We then 
sought to expand our analysis of their expression in GBM tis- 
sues through The Cancer Genome Atlas (TCGA), a publicly 
available repository which has accumulated comparative 
genomic hybridization, gene expression, and DNA methyla- 
tion analyses for 218 GBM samples [24]. Using this database, 
coexpression of MELK in GBM tumor samples was con- 
firmed (/j<.0001) with FOXMl expression. Immunocyto- 
chemistry with patient-derived GBM NS (GBM30) exhibited 
colocalization of FOXMl with MELK within individual cells 
(72.5% ± 2.4% colocalization in Fig. 3B). 

Since FOXMl and MELK are expressed in the same cells, 
we next investigated whether levels of MELK influence 
FOXMl expression and vice versa (Fig. 3C; Supporting 
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Figure 2. FoxMl expression is markedly elevated during gliomagenesis. (A): Representative images of immunohistochemistry for FoxMl (brown) 
in the SVZ of wild-type mice, SVZ of pretumorigenic Mut6 mice, and GBM-like tumors in cerebral cortex or normal areas in brain stem in Mut6 
mice. Original magnification is x 10. Scale bars = 100 /im. Magnified images are shown below. Original magnification is x40. Scale bars = 20 |im. 
(B): Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis for FoxMl expression in GBM tissues in Mut6 mice compared 
to normal side brain tissues derived from the same mice. Asterisks (*) indicate statistical significance by t test. Experiment repeated in three Mut-6 
mouse-derived GBM tissues. (C): Phase bright representative images of NS derived from either normal SVZ or GBM-like tumors developed in Mut6 
mice. Original magnification: x 10. Scale bars = 100 /im. Lower panel indicates qRT-PCR data for FoxMl expression in NS derived from GBM-like 
tumors (GBM NS) in Mut6 mice or nornial SVZ (normal NS). Asterisk (*) indicates statistical significance by / test. Experiments were done in tripli- 
cate and repeated three times independently. (D): Microarray analysis of relative FOX genes family expression in GBM (« = 30) compared to normal 
astrocytes cultures (n = 4). Expression of astrocytes was normalized as one. Expression of FOXMI is shown in red. (E): Relative expression of 
FOXMI of 10 GSCs samples compared to normal astrocytes cultures. FOXMI expression of normal astrocytes was normalized as one. (F): Relative 
mRNA expression of FOXMI and MELK in GBM157 and GBM528 CD133 (+) and CD133 (-) cells collected after cell sorting (n = 3). (G): Repre- 
sentative immunocytochemistry of GBM NS and GBM cells propagated in serum-containing medium (SPGCs) derived from GBM30. Cells were 
double-stained for FOXMI (red) in combination with one of the following neural progenitor cell markers Nestin, SOX2 and differentiation markers 
TuJl, GFAP (green). Hoechst dye is used for nuclear staining (blue). Original magnification: x40. Scale bars = 20 |im. Right panel indicates of the 
proportions of FOXMI positive cells coexpressing SOX2, Nestin, GFAP, and TuJl positive cells. For quantification, GBM30 sample were used and 
the experiment was repeated four times. (H): Flow cytometry analysis for the expression of FOXMI, SOX-2, and GFAP using GBM NSs and SPGCs, 
both of which are derived from GBM30. Abbreviations: FITC, fluorescein isothiocyanate; FSC, forward scatter; GFAP, glial fibrillary acidic protein; 
GBM, glioblastoma multiforme; NS, neurosphere; SPGC, serum-propagated GBM cells; SVZ, subventricular zone. 



Information Fig. S8). Both RT-PCR and Western blot demon- 
strated that MELK expression is not affected by FOXMI trans- 
genic expression or silencing. Likewise, RT-PCR analysis dem- 
onstrated that neither MELK overexpression nor knockdown 
affected FOXMI expression. We then asked whether FOXMI 
physically interacts with MELK. HEK 293 cells overexpressing 
both f /flg-tagged MELK and the fusion protein of FOXMI and 
EGFP or EGFP alone (control) were subjected to EGFP trap, 
followed by Western blot with MELK or Flag antibody (Fig. 
3D, upper panel; Supporting Information Fig. S9). Cells co- 
overexpressing MELK and FOXMI exhibited a stronger band 
for MELK than the control samples, implying a physical inter- 



action between MELK and FOXMI protein. This protein-pro- 
tein interaction is likely dependent on MELK kinase activity, as 
the band intensity for the co-immunoprecipitation (IP) of the 
catalytically dead MELK mutant (DI50A) was substantially 
diminished. Furthermore, we confirmed the direct interaction of 
these two proteins with an in vitro interaction study using 
recombinant purified FOXMI and the purified N-terminal 
MELK-fragment (amino acid 1-340) but not with the C-termi- 
nal fragment (not shown). 

We next asked whether recombinant MELK phosphorylates 
FOXMI (Fig. 3D, lower panels). Purified recombinant FOXMI 
was incubated with or without the purified kinase domain of 
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Figure 3. FOXMl is a substrate for MELK. (A): Comparison of MELK and FOXMl expression profile (Affymetrix Human Genome U133A 
Array) indicates statistically significant correlation of the expression of these two genes in grade III glioma (left) (n = 24) and in GBM (middle) 
(« = 56). The Cancer Genome Atlas (TCGA) analysis of MELK and FOXMl expression profile in 218 GBM patient samples indicates statistical 
significance (right panel). (B): Representative images of immunocytochemistry with GBM30 neurospheres for FOXMl (red), MELK (green). 
Hoechst dye for nuclear staining (blue). Original magnification: x40. Scale bar = 20 /(m. For quantification, GBM30 samples were used and the 
experiment was repeated four times. (C): Reverse transcription polymerase chain reaction (RT-PCR) analysis for MELK and FOXMl expression 
in GBM neurospheres treated with siRNA targeting three different sequences for FOXMl (left). RT-PCR analysis for MELK and FOXMl expres- 
sion in GBM neurospheres treated with different doses of siRNA targeting MELK (right) n = 5. (D): Upper panel: Overexpression of Empty- 
EGFP (control) + MELK-Flag (Flag-MELK), EGFP-FOXMl + MELK-F/ag, or EGFP-FOXMl + MELK D150A-f/a^; plasmids in HEK293 
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rylation of FOXMl by the kinase domain of MELK (1-340). Lane 1: FOXMl + ATP-Mg (no kinase); lanes 2 and 6: FOXMl + ATP-Mg + 
MELK (1-340); lane 3: FOXMl + MELK (1-340) (no ATP-Mg); lane 4: FOXMl + ATP-Mg + CDK2/CyclA; lane 5: FOXMl + ATP-Mg + 
CDK2/CyclA + roscovitine; lane 7: FOXM1+ ATP-Mg + MELK1-340(D150A). Lower panels: Coomassie staining of the samples subjected to 
autoradiography. Abbreviations: GBM, glioblastoma multiforme; IP, immunoprecipitation. 



MELK in an in vitro phosphorylation assay. In parallel, we used 
purihed CDK2/cyclin A kinase, a known FOXMl kinase, as a 
positive control [37]. The autoradiogram clearly shows in vitro 
phosphorylation of FOXMl by MELK and CDK2/cyclinA. In 
addition, FOXMl phosphorylation was not observed when incu- 
bated with kinase dead mutant MELK D150A. The same amount 
of FOXMl was present in all conditions, as is shown on the cor- 
responding Coomassie staining. Collectively, these data indicate 
that MELK directly interacts with FOXMl and thereby, phos- 
phorylates FOXMl in a kinase dependent manner. 



MELK Activates FOXMl Transcriptional Activity 
Leading to Upregulation of Mitotic Gene Expression 

To understand more in depth the molecular mechanism of 
MELK-dependent FOXMl signaling, we investigated whether 
MELK regulates FOXMl transcriptional activity using the 6x 
FOXMi-TATA-luciferase reporter plasmid. Transgene expres- 
sion of wild-type FOXMl alone activated its transcriptional 
activity, suggesting self-activation in GBM30 GSCs (Fig. 
4A). Transgenic expression of wild-type MELK significantly 
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Figure 4. MELK phosphorylates FOXMl and regulates FOXMl activity, leading to upregulation of mitotic gene expression. (A): Graph indi- 
cating the FOXMl promoter activity in 293T cells transfected with the 6x FOXMl TATA-luciferase plasmid together with expression vectors 
encoding wild-type FOXMl and increasing amounts of plasmids encoding either WT MELK or kinase-dead mutant form of MELK (D150A). 
The experiment was performed in triplicate in 96 well plates and repeated three times independently. (B): Relative mRNA expression levels of 
Survirin, CyclinBl, CDC25B, and Aurora B by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in GBM30 cells trans- 
fected with OFF, FOXMl, FOXMl + WT MELK, and FOXMl + MELK D150A mutant. qRT-PCR was performed in triplicate and repeated 
three times independently. Abbreviations: EGFP, enhanced green fluorescent protein; WT, wild type. 



enhanced FOXMl self-activation in a concentration-dependent 
manner (Fig. 4A). On the other hand, transgenic expression of 
the kinase dead mutant MELK D150A eliminated FOXMl 
self-activation in a concentration-dependent manner. Since 
FOXMl is known to regulate the transcriptional network of 



genes essential for mitotic progression (e.g., Aurora B kinase, 
CENPA, Siirvivin, and CyclinBl), we examined whether the 
FOXMl -MELK complex had any effect on mitotic regulators 
in GBM30 GSCs. Cotransfection of FOXMl with MELK 
resulted in a marked elevation of the key mitosis genes (Fig. 
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4B) [38]. In contrast, MELK D150A did not influence 
FOXMl -mediated mitotic transcriptional expression. Taken 
together, MELK kinase activity regulates phosphorylation and 
self-activation of FOXMl, leading to enhanced mitotic gene 
expression in GSCs. 



MELK-Driven Phosphorylation of FOXMl Is PLKl 
Dependant 

We then sought to determine the mechanisms that regulate 
FOXMl activation by MELK in GSCs. A recent study by Fu 
et al. [39] demonstrated that, similar to MELK, PLKl binds 
and phosphorylates FOXMl leading to FOXMl activation and 
increased gene expression required for mitotic progression in 
hepatocellular carcinoma cells. We therefore sought to dissoci- 
ate the molecular signaling mechanisms of MELK-driven 
FOXMl phosphorylation and PLKl -driven FOXMl phospho- 
rylation in GSCs. First, we used the PLKl specific inhibitor 
ON01910 together with overexpression of FOXMl, MELK, 
and/or PLKl and evaluated FOXMl promoter activity (Fig. 
5A). After confirming the elimination of PLKl -driven FOXMl 
activation by ONO1910, we found that ONO1910 also elimi- 
nates FOXMl activation driven by MELK in a dose-dependent 
manner (Fig. 5A). In order to exclude the possibility of an off- 
target effect by pharmacological inhibition, we then tested the 
effect of the PLKl dominant negative mutant (PLKl DN) on 
MELK-mediated FOXMl activation. Similar to the results with 
the pharmacological inhibition, forced expression of PLKl DN 
eliminated MELK-driven FOXMl activation in GBM30 GSCs 
(Fig. 5B). To further confirm these results, we used the consti- 
tutive active form of FOXMl (FOXMl EE as PLKl -driven 
phosphomimics) to assess whether MELK D150 mutant is able 
to eliminate FOXMl EE-driven FOXMl activation in 293T and 
GBM30 GSCs (Fig. 5C). Unlike the MELK D150A-induced 
elimination of FOXMl wild-type-driven self-activation, MELK 
D150A failed to diminish FOXMl activation driven by 
FOXMl EE in both 293T cells and GBM30 GSCs. Taken to- 
gether, these data suggest that MELK-regulated FOXMl activa- 
tion is dependent on PLKl in GSCs. 

We next asked whether MELK-driven FOXMl activation is 
regulated by either the initial binding event of PLKl to 
FOXMl or the secondary trans-phosphorylation of FOXMl fol- 
lowing PLKl-FOXMl binding. Specifically, Ser 715 and Ser 
724 of the FOXMl protein were mapped to be the major trans- 
phosphorylation sites and Thr-596 and Ser-678 were mapped to 
be the PLKl binding sites [39]. A single mutation of Ser 715 or 
Ser 724 partially reduces the FOXMl self-activation, whereas a 
transgenic expression of the double mutation results in com- 
plete abolishment of FOXMl trans-activation [39]. Mutation of 
one of the two residues resulted in only a partial decrease in 
MELK-enhanced FOXMl self-activation in both 293T cells 
and GBM30 GSCs (Fig. 5D). On the other hand, introducing 
the double mutation led to a major reduction in MELK-depend- 
ent FOXMl transcriptional activity nearly to the level resulting 
from the loss of MELK kinase activity induced by the D150A 
mutation (Fig. 5D). We then performed the same experiments 
with the FOXMl mutant vectors lacking the ability of PLKl to 
bind FOXMl either partially (single mutation vectors) or com- 
pletely (double mutation vector) (Fig. 5E). Transgenic expres- 
sion of the double mutation vector (FOXMl 596/678), but not 
the single mutation vectors (FOXMl T596A and S678A), 
almost completely abrogated MELK-driven FOXMl activation 
both in 293T cells and GSCs. Taken together, these data indi- 
cate that the MELK-driven FOXMl activation is dependent on 
both the initial binding and subsequent trans-phosphorylation 
events by PLKl in GSCs. 



Siomycin A Abrogates FOXMl-MELK Interaction 

Previously, the thiazole antibiotic Siomycin A was identified 
as a FOXMl inhibitor by diminishing its protein and mRNA 
abundance, and we recently discovered that Siomycin A also 
reduced MELK expression in GSCs in vitro [9, 18, 40]. Con- 
sistent with these findings, flow cytometry demonstrated that 
Siomycin A reduced the proportions of FOXMl (+) and 
MELK (+) cells in GBM NS in a dose-dependent manner 
(Fig. 6A). We then investigated the effect of Siomycin A on 
FOXMl activity driven by FOXMl and/or the FOXMl-MELK 
complex (Fig. 6B). We treated GBM30 GSCs with 0.1 /(M or 
0.5 ^M of Siomycin A after transfection. Similar to the protein 
expression, FOXMl activity was almost abolished, declining 
to basal levels at 48 hours post-treatment with 0.5 /(M of Sio- 
mycin A, suggesting that Siomycin A abrogates FOXMl activ- 
ity driven by FOXMl itself as well as MELK-driven FOXMl 
phosphorylation. Next we evaluated whether restoring 
FOXMl and MELK rescues Siomycin A-induced mitotic 
arrest in GSCs (Fig. 6C). In agreement with previous studies 
using other cancer cells, GSCs increased the proportion of 
cells in the G2/M phase when treated with 0.5 pM of Siomy- 
cin A (cefls in G2/M phase: 18.72%-34.23%) [40, 41]. This 
data indicated that Siomycin A induces G2/M arrest of GSCs. 
Overexpression of FOXMl with MELK, but not FOXMl or 
MELK alone, rescued the Siomycin A-induced G2/M arrest of 
GBM30 GSCs (cells in G2/M phase: 34.23%-21.74%). These 
data provided supporting evidence for the critical role of the 
FOXMl-MELK complex in mitotic progression of GSCs. 

Since FOXMl levels were markedly higher in GSCs com- 
pare to NPCs (Fig. 2), we reasoned that the dependence of 
GSCs on FOXMl would be greater than that of NPCs. 
Indeed, when we treated mouse and human GSC cultures 
(derived from Mut 6 mice and GBM30, respectively) and 
NPC cultures (derived from wild-type mice and a 16- week- 
old fetus-16wf, respectively) with varying doses of Siomycin 
A, GSCs derived from both species were 15-fold more sensi- 
tive to Siomycin A than NPCs, suggesting that pharmacologic 
inhibition of FOXMl signaling by Siomycin A has stronger 
impact on in vitro growth of GSCs than NPCs (Fig. 6D). 

Combined Treatment of Temozolomide with Siomy- 
cin A on GSCs-Derived Mouse Tumors Yields Better 
Survival than Monotherapy with Temozolomide 

Accumulated evidence suggests that GSCs are relatively re- 
sistant to Temozolomide chemotherapy partly due to elevated 
expression of the drug resistance genes compared to non- 
GSCs [42^5]. We found that Temozolomide treatment of 
GBM30 GSCs paradoxically increases the expression of 
FOXMl and MELK in a dose-dependent manner (Fig. 7 A). 
These findings are consistent with the hypotheses that Temo- 
zolomide treatment spares FOXMl and MELK expressing 
GSCs or promotes the appearance of a GSC-like phenotype. 

We then investigated the effect of combination treatment 
with Temozolomide and Siomycin A on GSCs-derived mouse 
tumors in vivo (Fig. 7B-7E). Xenografting GBM30 NS 
resulted in the formation of GBM-like tumors in immunocom- 
promised mouse brains. Similar to human GBM, the 
xenograft-derived tumors exhibited massive intratumoral hem- 
orrhage, necrosis, and densely packed tumor cells with hyper- 
chromatic nuclei surrounding tumor vessels (Fig. 7B). 
Immunohistochemistry of mouse tumors at 2 days post- 
Siomycin A injection demonstrated a significantly reduced 
fraction of Ki-67(+) proliferating cells compared to control- 
treated mouse tumors (Fig. 7C). In order to determine the 
effect of Siomycin A treatment on GSCs in vivo, we cultured 
NS forming cells of treated and control tumors. Tumors 
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Figure 6. SM abrogates MELK-driven FOXMl activity and inhibits growth of glioma stem-like cells (GSCs) but not normal progenitors in mouse 
and human. (A): Flow cytometry analysis on GBM30 spheres treated with either DMSO or SM. The upper panels indicate the proportions of 
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cells in each condition. Experiment was repeated three times and similar expression pattern was observed. (B): Graph indicating the FOXMl pro- 
moter activity in indicated cells transfected with the plasmids encoding for FOXMl wild-type (WT) with MELK WT or kinase dead D150A. The 
cells were then treated with indicated doses of SM after transfection for 48 hours, followed by detection of FOXMl activity by luciferase assay. 
(C): Graph showing cell cycle analysis with flow cytometry of GBM30 GSCs. The cells were treated with either control (DMSO 1%), SM (500 
nM), FOXMl wild-type overexpression together with SM (500 nM), MELK wild-type overexpression together with SM (500 nM) or MELK, and 
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DMSO, dimethyl sulfoxide; FITC, fluorescein isothiocyanate; FSC, forward scatter; GBM, glioblastoma multiforme; SM, Siomycin A. 



assayed at 2 days post-Siomycin A injection had significantly 
fewer NS-forming cells compared to the control tumors (Fig. 
7D). These data suggest that Siomycin A has a potent inhibi- 
tory effect on survival and proliferation of GSCs in vivo. Sub- 
sequently, we assessed if Siomycin A treatment has any sur- 
vival benefit on Temozolomide-treated mouse tumors. When 
the tumors were treated with Temozolomide alone, median 
survival of tumor-bearing mice was prolonged from 11.5 days 
to 22.5 days (Fig. 7E). The combined treatment with Temozo- 
lomide and Siomycin A exhibited further benefit on mouse 
survival (median survival to 28 days) (DMSO control vs. 
Temozolomide+Siomycin A: p< .001). 

Finally, we used the tumor slice culture method to assess 
Siomycin A treatment on surgical GBM specimens [46]. 
Patient-derived GBM slice cultures were treated with DMSO, 
Temozolomide, or Siomycin A, and the effects were measured 
16 hours post-treatment. H&E staining of the slice cultures 
demonstrated that the procedures did not destroy the cytoarchi- 



tecture of GBM tumors (Supporting Information Fig. SIO). 
Immunohistochemistry with a proliferation marker, Ki67 dem- 
onstrated vast numbers of proliferating tumor cells in DMSO- 
treated samples, but not in Temozolomide- or Siomycin A- 
treated tissues (Supporting Information Fig. Sll). In turn, 
immunohistochemistry for apoptosis marker activated Caspase- 
3 displayed significantly higher number of apoptotic cells in 
Siomycin A-treated samples, but not in DMSO- or Temozolo- 
mide-treated cultures (Supporting Information Fig. Sll). 



Discussion 



Molecular signaling between protein kinases and transcription 
factors plays vital roles in tumor development and mainte- 
nance [47]. Here, we demonstrate the first evidence that the 
transcription factor/oncogene FOXMl forms a protein 
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complex with a serine/threonine kinase MELK and that 
FOXMl serves as a substrate of MELK in cancer cells, and 
MELK-regulated FOXMl phosphorylation controls FOXMl 
activity and induces the expression of downstream mitotic 
regulators. The transcription factor FOXMl is a master regu- 
lator for cell cycle progression and is overexpressed in a num- 
ber of human cancers including GBM [48]. The protein kinase 
MELK is also abundantly expressed in various cancers includ- 
ing GBM and plays a pivotal role in survival of cancer cells 



and cancer stem cells [12], Recent studies have shed light on 
FOXMl signaling in GBM and GSCs [10]. FOXMl interacts 
with the proiTioter of the VEGF gene regulating its activation, 
which contributes to GBM tumor angiogenesis [49]. Another 
study by Zhang et al. elegantly showed that FOXMl is a 
downstream component of canonical Wnt-signaling and 
directly binds to /?-catenin, inducing its nuclear translocation 
and transcriptional activation [39]. This study demonstrates 
the first evidence that FOXMl directly interacts and is 
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phosphorylated by MELK in GSCs, further providing evi- 
dence for a role of FOXMl signaling in GBM. Our data indi- 
cate that both phosphorylation of FOXMl protein and tran- 
scriptional activity of FOXMl depend on MELK kinase 
activity. Intriguingly, FOXMl activation driven by MELK is 
dependent on both the priming event of the complex forma- 
tion of FOXMl with PLKl and the subsequent trans-activa- 
tion of FOXMl by PLKl. 

Despite the advances presented here, open questions still 
remain. First, it is not clear whether the MELK/FOXMl pro- 
tein complex plays a positive role in cancers or cancer stem 
cells in other organs. Second, it is undetermined whether 
MELK-mediated signaling is associated with the other pro- 
teins that regulate, or are regulated by, FOXMl in cancer 
cells. Future studies are needed to elucidate these questions. 

The cell of origin for GBM is still debatable, although the 
concept has been evolving. Here, we used a mouse model of 
gliomagenesis by ablation of key tumor suppressor genes in 
NPCs causing formation of GBM-like tumors in cerebral cor- 
tices. Our data indicated that FOXMl signaling is operative 
in NPCs of the developing brains in vivo. Immunohistochem- 
istry exhibited that abundance of FOXMl protein is progres- 
sively elevated as cells progress from the neural progenitor 
stage through pretumorigenic progenitors to GSC. Sensitivity 
of GSCs to Siomycin A was markedly higher than that of 
NPCs. Taken together, these data indicate that GSCs are more 
dependent on FOXMl signaling compared to NPCs. However, 
it remains unknown whether activation of FOXMl signaling 
alone is sufficient for gliomagenesis and whether FOXMl sig- 
naling plays a critical role in oligodendrocyte precursor- 
derived gliomagenesis. Future studies could help address these 
questions. 

An important therapeutic implication of the present data is 
that combined therapy of mouse xenografted tumors with 
Temozolomide and Siomycin A resulted in diminished NS- 
forming GSCs in tumors and provided an additional benefit 
on mouse survival. Although eradication of cancer stem cells 
appears to be essential for the cure of cancers, recent studies 
also suggest that non-cancer stem cells acquire cancer stem 
cell phenotype when challenged by stressors such as ionizing 
radiation and chemotherapies [50]. To achieve cure of a can- 
cer, eradication of the rare population of tumor cells that are 
considered as existing cancer stem cells may not be sufficient; 
instead, combination of non-cancer stem cell-targeted thera- 
pies and cancer stem cell-targeted therapies appears to be 
mandatory. Our data indicate that Temozolomide treatment 
increased, while Siomycin A decreased, both FOXMl and 
MELK expression in GSCs. One interpretation of these data 
is the preferential eradication of non-stem tumor cells and 
subsequent enrichment of GSCs after treatment (selection of 
therapy-resistant tumor cell population; the clonal evolution 
theory). Alternatively, Temozolomide may induce phenotypic 



changes of the treated GBM cells and may increase FOXMl 
and MELK expression (the molecular evolution theory). It is 
also possible that both theories are true. Future studies will 
address this open question. 



Conclusion 



In this study, we demonstrated that the oncoprotein FOXMl 
is regulated by the mitotic kinase MELK for its phosphoryla- 
tion and autoactivation via a direct interaction in GSCs. 
Although FOXMl is restrictedly expressed in NPCs in the 
normal brain and plays a role in NPC growth, a substantial 
elevation of FOXMl during gliomagenesis and elevated sensi- 
tivity of GSCs over NPCs to pharmacological inhibition of 
FOXMl signaling indicate that GSCs are more dependent on 
FOXMl signaling for their survival and growth. Thus, target- 
ing the ability of FOXMl to form a protein complex with 
MELK may represent a potential therapeutic benefit for 
GBM. The results presented here may help to gain further 
insights into the biology of GSCs as well as molecular mech- 
anism of tumorigenesis and therapy resistance in GBM. 
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